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SECTION  I 
INTRODUCTION 


The  Eglin  Aeroballistics  range  is  an  enclosed  structure  for  recording 
projectile  flight  paths,  Reference  1.  The  projectiles  are  in  various  shapes  and 
may  travel  at  speeds  up  to  MACH  6.  They  are  imaged  at  regular  intervals  along 
the  flight  path  so  that  an  aerodynamic  analysis  can  be  conducted.  At  each 
recording  station,  orthogonal  image  pairs  are  collected.  Each  station  is  out¬ 
fitted  with  two  4  by  5  cameras  as  shown  in  Figure  1.  Various  support  electronics 
accompanying  the  camera  provide  shutter  control  and  time-recording  functions. 

The  projectile  is  illuminated  by  a  Hi-Voltage  model  SS55P  spark-gap  flash 
positioned  adjacent  to  each  camera.  The  flash  duration  is  approximately 
0.3  ps.  The  flash  can  be  considered  a  point  source  of  light;  therefore,  the 
projectile  casts  a  well-defined  shadow  on  the  special  reflective  surface  (3M 
No.  7610)  on  the  opposite  wall.  The  flash  is  activated  at  a  predetermined  time 
after  the  projectile  passes  an  infrared  linear  array  positioned  near  the  camera. 

Our  task  is  to  investigate  replacing  the  existing  film  cameras  with 
electronic  imaging  (El)  devices  to  reduce  the  delays  and  the  labor-intensive 
operation  associated  with  the  present  facility.  The  conceptual  design  requires 
that  the  resolution  of  the  existing  system  be  maintained  and  the  cameras  be 
replaced  without  perturbing  the  rest  of  the  system. 


Figure  1.  Typical  Shadowgraph  Station 


SECTION  II 

ANALYSIS  OF  FACILITY 


Physical  and  temporal  restrictions  imposed  by  modern  aeroballistics  have 
pushed  imaging  technology  to  the  point  where  special  photoconductive  surfaces 
and  high-speed  support  electronics  are  dictated.  Specifications  for  these 
devices  can  be  formulated  by  a  methodical  analysis  of  critical  parameters  and 
their  interaction.  In  terms  of  system  theory,  system  transfer  functions  and 
state  equations  can  be  used  in  optimal  coupling  of  devices  to  maximize  system 
performance.  Application  of  these  methods  to  El  at  Eglin  is  described  in  this 
report . 

The  El  system  must  autonomously  collect  and  permanently  archive  an  image. 
It  must  provide  capabilities  equivalent  to  existing  hardware,  Reference  2; 
performance  has  been  defined  relative  to  the  existing  equipment.  The  El 
system  must  be  capable  of  acquiring  an  image  in  0.3  ps.  An  air-gap  flash 
will  illuminate  the  object  field  with  sensor  spectral  optimization  provided 
for  low  light  requirements.  The  system  resolution  must  be  consistent  with  the 
existing  equipment. 

From  the  onset,  it  was  clear  that  a  spatial,  temporal,  and  spectral 
characterization  of  the  system  was  needed  to  accurately  specify  El  hardware; 
more  fundamentally,  such  a  characterization  was  needed  to  determine  if  the 
project  was  realizable  with  current  state-of-the-art  devices.  Therefore, 
theoretical  analysis  of  the  system  was  done  to  justify  the  underlying  princi¬ 
ples  of  the  proposed  El  system.  The  analysis  resulted  in  a  favorable  outlook 
for  the  proposed  system.  The  next  step  was  to  analyze  the  system  by  means  of 
empirical  data.  This  analysis  consisted  of  two  primary  determinations:  the 
system  transfer  function;  and,  the  spectral  output  of  the  light  source. 

First,  system  resolution  was  defined  in  terms  of  signal-to-noise  (S/N) 
ratio,  as  compatible  with  commercial  sensor  data.  Second,  spectral  response 
of  the  system  was  quantified,  which  was  necessary  for  optimizing  the  spectral 
sensitivity  of  the  El  sensor  to  the  existing  light  source.  These  two  steps 
are  covered  in  this  section. 

Theoretical  analysis  established  initial  problems  and  concepts  in  El. 
General  digital  El  resolution  limitations  were  evaluated  along  with  image 
processing  techniques  for  a  posteriori  error  correction.  In  addition,  the 
study  in  Sec.  II. 6.  established  inherent  limitations  in  both  Kodak  Royal  X-Pan 
film  and  El  devices. 
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1 .  IMAGE  ACQUISITION 

The  images  to  be  analyzed  are  two-dimensional  projections  of  fast 
changing  three-dimensional  systems.  Therefore,  we  needed  an  accurate  time 
history  of  a  projectile  in  flight  in  three  dimensions.  Recording  orthogonal 
projections  at  various  sample  times  during  the  flight  provided  this  accurate 
history.  Complete  three-dimensional  information  pertaining  to  the  attitude  of 
the  projectile  can  be  obtained  if  the  projections  are  orthogonal  to  each  other 
and  if  the  imaging  system  has  sufficient  resolution. 

A  digital  image  is  a  discrete  representation  of  a  continuous  real  image. 
Enough  samples  of  a  continuous  image  must  be  taken  to  preserve  the  fine  detail 
Cor  resolution) .  Field-of-view  and  resolution  concepts  must  be  examined  to 
determine  minimal  specifications  of  the  imaging  system.  A  field  of  view  is 

normally  an  image  that  consists  of  an  array  of  numerical  points  (N  x  N  ) . 

P  P 

The  field  of  view  must  be  large  enough  to  contain  the  projectile  each  time  it 

is  fired,  since  it  is  uncertain  where  the  projectile  will  be  in  the  field  of 

view.  Enlarging  the  field  of  view  without  adding  more  points  (that  is,  making 

N  larger)  lessens  the  detail  that  can  be  measured  in  the  vicinity  of  the 
P 

projectile;  therefore,  it  is  desirable  that  the  projectile  image  fill  up  the 
field  of  view  as  much  as  possible. 

With  the  field-of-view  and  resolution  relationships  known,  the  amount  of 
fine  detail  in  the  projectile  image  to  be  preserved  must  be  specified.  Using 
the  two-dimensional  Nyquist  sampling  criteria  that  require  sampling  at  spatial 
frequency  twice  that  of  the  finest  detail  satisfies  this  requirement.  If 
spatial  sampling  frequency  is  not  high  enough,  two-dimensional  aliasing  will 
produce  artifacts  and  moire  effects  in-  a  reconstructed  image. 

Only  careful  examination  of  all  parameters  in  an  imaging  system  and  of 
the  results  desired  can  dictate  the  actual  specifications  to  be  built  into  any 
imaging  system.  However,  certain  crude  estimates  of  imaging  system  perform¬ 
ance  can  be  made  for  a  set  of  assumptions  that  pertain  to  the  system.  Con¬ 
sider  a  projectile  image  whose  long  axis  fills  up  fraction  f  of  the  field  of 
view  of  an  Np  x  N^  image.  Further,  assume  that  the  projectile  image  is 
sampled  at  the  two-dimensional  Nyquist  spatial  frequency.  Then  geometric 
considerations  alone,  given  uncertainties  in  the  location  of  the  end  points  of 
the  long  a.cis,  will  introduce  a  minimum  error  associated  with  any  measure¬ 
ment.  Any  angle,  measured  for  the  long  axis,  can  be  in  error  by  +0,  where 


Other  errors  in  measurement,  such  as  improper  Nyquist  rate,  will  add  to  this 
minimum  error.  The  relationship  between  image  accuracy  and  variation  of  N^ 
and  f  is  clearly  shown  in  Table  I. 


f 

TABLE  1.  LONG-AXIS  ANGLE  MEASUREMENT  ERROR 
(3  degrees) 

N 

_ 

1000 

500 

250 

100 

0.05 

2.4 

5.0 

11.0 

34.0 

0.10 

1.2 

2.4 

5.0 

14.0 

0.25 

0.46 

0.93 

1.9 

5.0 

0.50 

0.23 

0.46 

0.93 

2.4 

A  POSTERIORI 

IMAGE  PROCESSING 

Image  processing  is  the  manipulation,  usually  by  digital  computer,  of 
the  acquired  image  with  intent  to  detect  or  enhance  some  feature.  For  a 
projectile,  image  processing  should  dynamically  measure  the  attitude  of  the 
model  in  two  orthogonal  projection  planes.  Various  image  processing  opera¬ 
tions,  such  as  edge  detection  and  smoothing  to  eliminate  moire  effects,  will 
be  required. 

a.  Correcting  System  Distortion  and  Noise.  Camera  lenses  commonly 
exhibit  nonlinear  distortion  both  in  the  radial  distance  from  the  lens  mid¬ 
point  axis  and  in  the  distance  of  the  object  from  the  image  plane.  Figure  2. 

In  addition,  the  detector  can  distort  the  image  by  nature  of  the 
physical  recording  process  and  add  noise  because  of  digitization. 
Characterizing  the  lens  and  the  detector  mathematically  can  remove  the 
distortion  effects  to  significantly  improve  image  resolution. 

b.  Matnematical  Model.  Bach  component  of  the  camera  system  affects  the 
resulting  image  and  can  be  treated  separately.  The  camera  lens  and  the 
television  sensor  or  array  detector  act  on  the  image-forming  light. 

Figure  3.  In  addition  to  the  effect  of  the  detector  on  the  image,  the 
sensors  s{b}  may  have  a  signal-dependent  noise  effect  n^  on  the  image. 
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OBJECT 


Figure  2.  Image  Geometry  Model 


n3(x,y) 

SIGNAL-DEPENDENT  NOISE 


f 

OBJECT  PLANE 
RADIANT 

ENERGIES 


h 


IMAGE 

FORMATION 

DETECTOR 

(LENS) 

b  (x 

-y) 

s(b) 

IMAGE 

PLANE 

RADIANT  ENERGIES 


r(x,y) 
DETECTOR 
RESPONSE 


RECORDED 
IMAGE 

PLUS  NOISE 
g  (x  ,y )  =  r(  x  ,y )+  n2  (x^ ) 

+n3<xfy) 


n2  (x,y) 

NOISE  PROCESS 


Figure  3.  Transfer  Function  Model  of  the  Imaging  Process 


Quantization  and  electrical  noises  are  added  noise  sources  n^  that  also  affect 
the  image . 

c.  Image  Reconstruction.  The  recorded  imag>  from  a  camera  is  an  array  of 

2 

H  data  entries.  Using  the  lens  point-spread-function  matrix,*  we  can 
represent  the  lens  output  in  state  space  as 

g  -  IHlf  . 

where  i  is  the  image  plane  matrix  with  I  by  H  elements  contiguously  arranged 
into  a  vector  and  f  is  the  H  by  II  matrix  digital  representation  of  the 
object  plane  contiguously  arranged  into  a  vector.  H  is  the  frequency  domain 
transfer  function,  and  its  dimensions  are  the  square  of  the  number  of  ele¬ 
ments  in  the  image  plane  (in  practice  many  of  the  elements  will  be  zero  or 
contain  no  information) .  The  transfer  function-in-state  space  of  the  entire 
system  (Figure  3)  is 

g  -  s{ [H]f }  +  d{ [H]  f  }  +  n2  . 

How  the  image  f '  can  be  restored  by  the  nonlinear  operator  T. 

?'  =  T [£]  -  T{s([H]f)  +  d( [H]f )  +  52J  . 

Fortunately,  television  cameras  are  equipped  with  corrector  circuits  to  cancel 
signal-dependent  noise  produced  by  the  sensors;  therefore, 

f*  =  T{s([H]f)  +  n2J  . 

Then  the  natural  choice  for  the  restoration  process  is 
T  =  [Hl^s-1  ; 

*The  lens  point-spread-function  matrix  can  be  obtained  at  Los  Alamos 
facilities  or  from  commercial  optics  laboratories. 


therefore. 


f'  *  iHrV^sdHlf)  +  n2)  ,  and 

f*  »  f  +  tHl”1s”152  . 

The  result  is  a  restored  inage  f '  exactly  like  the  original  object  plus  the 
noise  acted  on  by  the  inverse  detector  and  lens  point-spread  functions.  The 
noise  is  nominally  small  thereby  reducing  the  significance  of  this  term. 

The  camera  lens  may  significantly  distort  the  object  relative  to  the 
distance  from  the  image  plane.  We  can  use  the  above  equations  to  correct  this 
distortion,  but  we  must  know  the  projectile  size  and  object-to-image  scale 
factor  to  calculate  the  distance  from  the  actual  image.  The  scale  factor 
would  be  a  constant  determined  when  the  camera  was  installed.  Projectile  size 
would  be  an  input  parameter  fed  to  the  computer,  and  the  image  plane  size 
would  be  calculated  by  the  computer. 

d.  Projectile  Position.  The  exact  position  of  the  projectile  tip  is 
difficult  to  determine  directly  because  of  the  finite  number  of  points  that 
make  up  the  projectile  image  when  the  dominant  spectral  content  of  the  shape 
of  the  projectile  tip  exceeds  the  Hyquist  frequency  limit.  The  Vyquist 
frequency  limit  can  be  calculated  from  the  general  expression  of  uniform  grid 
sampling,  and  then  an  expression  can  be  derived  for  the  minimum  radius  of 
curvature.  Hote  that  the  dimension  units  for  frequency  are  radians  per 
distance  using  the  grid  sampling  technique  (Figure  4): 


g  (mAx.ndy)  *  X  X  g(x,y)  4(x  -  mAx.y  -  nAy)  , 
s 

m  n 

where  4  is  the  Dirac  delta  function.  The  Fourier  transform  of  the  image  is 
(u  *  radians  per  distance,  i  *  -1) 


G(w  ,u  ) 
x  y 


/  fz (x.y)  exp  (-2wi((i>xx  +  w  y)l  dx  dy 


G  (u  ,U  )  a  G(u  ,W_)  *  S(w  .«  ) 


The  Nyquist  frequency  limit  is  then 


2  *  . 
w  ■  _  and 

xc  2  Ax 


2  « 

u  *  _  .  or 

yc  2  Ay 


u 


xc  Ax 


and  u 


yc  Ay 


When  we  approximate  the  radius  of  curvature  as  one-fourth  cycle, 

2*  _  2*  Ax 

4<a  “  4v/Ax  =  2 

xc 


Minimum  radius  of  curvature  is  then 


‘xmin 


Ax  . 

~Z 


ymin 
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For  projectile  tips  with  a  radius  of  curvature  less  than  r^^,  the  tip 
position  must  be  approximated.  The  tip  could  be  generated  from  statistically 
selected  points  on  the  projectile  perimeter  using  standard  predictive  estima¬ 
tion  routines. 

3.  SYSTEM  TRANSFER  FUNCTION 

The  system  transfer  function  was  acquired  by  using  empirical  data.  In 
this  way  the  effects  of  the  lens  system,  film  medium,  developing  process,  and 
illumination  source  are  encompassed  in  a  single  format  from  which  a  defini¬ 
tive  answer  to  system  resolution  can  be  drawn. 

The  basic  data  set  required  is  the  time  domain  impulse  response  f(t)  of 
the  system.  From  this  the  Fourier  transform.  Equation  (1)  can  be  applied  to 
get  the  transfer  function  F(x);  that  is, 

00 

F(x)  »  ?(f(t)J  »  /  f (t)exp[-j2wtx]dt  .  (1) 
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This  method  requires  use  of  an  impulse  as  the  system  input.  Because  this  is  a 
physically  unrealizable  entity,  a  step  function  was  used.  The  step  was 
produced  on  precision-ground  glass  using  special  in-house  photograph-masking 
equipment.  A  step  function  S(t)  can  be  related  to  an  impulse  function  4(t) 
by  Equation  (2). 


*S(t) 

dt 


Once  the  Fourier  transform  of  the  step  response  has  been  calculated,  the 
transfer  function  may  be  derived  by  relating  Equation  (2)  to  Fourier  space. 
Given  that 


f (x)  — -  F(t)  , 


the  inverse  Fourier  transform  is 


f(t)  -  f  F(x)exp( j2«xt)dx  . 


Differentiating  both  sides  gives 


df  (t) 
dt 


d 

—  f  F(x)exp(j2*xt)dx 

ac  -oo 
00 

/  (j2*x)  F(x)exp(j2*xt)dx  , 


and  comparing  Equation  (4)  with  Equation  (3)  clearly  shows  that  [df(t)]/dt  and 
(j2wx)  F(x)  are  a  Fourier  transform  pair.  Thus,  the  transfer  function  is 
obtained  by  dividing  the  Fourier  transform  of  the  step  response  by  j2*x. 

a.  Experimental  Considerations.  The  step  response  data  were  acquired 
using  existing  equipment  and  the  geometry  currently  used  for  high-speed  aero- 
ballistics  photography.  The  experiment  was  conducted  at  the  center  and  at  the 
four  extremities  of  the  field  of  view.  Data  from  the  four  off-axis  locations 
will  be  used  to  evaluate  lens  anomalies — a  factor  that  increases  as  the  radial 


Digitization  of  Step  Response  as  Recorded  on  Film  (Sampling 
interval  is  4  pa) 


•SCAN  DISTANCE  (>im) 

Figure  6.  Average  Response  to  Step  Input 


11 


distance  from  the  optical  axis  increases.  Our  experiment  was  conducted  at 
several  different  photographic  stations  with  multiple  exposures  at  each 
station  for  statistical  validity. 

Physically,  the  step  response  was  recorded  on  film  using  standard  equip¬ 
ment  and  operating  procedures.  The  data  were  prepared  for  computer  processing 
by  digitizing  the  photograph  using  a  Perkin-Elmer  scanning  microdensitometer. 
The  step  response  was  scanned  a  number  of  times  at  4-ym  increments 
(Figure  5);  then  an  average  of  all  scans  (Figure  6)  was  used  for  further 
processing. 

b.  Analysis.  A  plot  of  the  off-axis  system  transfer  function  [or  modu¬ 
lation  transfer  function  (MTF) ]  is  shown  in  Figure  7.  This  plot  indicates  the 
normalized  system  response  versus  frequency  in  cycles  per  millimeter  at  the 
object  field.  The  goal  here  is  to  provide  a  reasonable  estimate  of  the 
dimensions  of  the  smallest  object  that  can  be  resolved.  Studies  (Reference  3, 
page  119)  have  indicated  that  the  eye  normally  requires  an  image  contrast  of 
about  3  percent  to  resolve  light-to-dark  transitions.  On  a  normalized  MTF, 
this  contrast  value  is  the  spatial  frequency  at  which  the  response  has  fallen 
to  3  percent  of  the  passband  amplitude;  this  frequency  is  called  the  limiting 
resolution. 


FREQUENCY  (CYCLES/mm) 

Figure  7.  System  Transfer  Function  at  the  Extreme  Boundary  of  the  Field  of  View 


FREQUENCY  (CYCLES/mm) 


Figure  8.  System  Transfer  Function  at  the  Center  of  the  Field  of  View  (on  the 
optical  axis) 

The  center  field  frequency  response  is  shown  in  Figure  8.  Taking  the 
lower  frequency  limit  as  the  lowest  noise-free  point  yields  a  value  of  0.742 
cycles  per  millimeter.  Using  this  value  as  the  zero  reference  boundary,  we 
calculated  the  limiting  resolution  to  be  0.72  cycles  per  millimeter.  One 
industry  standard  for  system  specification  is  the  contrast  transfer  function 
(CTF)  that  has  square  wave  basis  functions  instead  of  sine  waves.  CTF  can  be 
calculated  from  the  HTF  by 

ctp  .  C(»>  -  4  [non  .  Mail ,  siai .  Sim  +  H£|ai ...]  . 

The  equation  yields  0.84  line  pairs  per  millimeter  for  C(H)  at  the  center  of 
the  field.  In  comparison,  the  comers  will  resolve  about  0.5  line  pairs  per 
millimeter.  This  difference  is  due  to  nonlinear  lens  effects  that  increase 
with  radial  distance  from  the  optical  axis.  Graphically,  this  difference  can 
be  seen  in  the  differing  cutoff  frequencies  between  the  off-axis  transfer 
function  of  Figure  7  and  the  on-axis  transfer  function  of  Figure  8. 


A  quick  way  of  verifying  the  CTF  resolution  derived  here  from  experimental 

data  is  to  calculate  the  theoretical  system  resolution  (II  ) .  This  is  deter- 

s 

mined  from  all  elements  in  the  optical  chain,  up  to  and  including  the  film. 


1  _L_ 

„  2  *  Z,  „  2  . 


i  Mi 


As  an  estimate,  the  major  components  can  be  used  in  this  evaluation. 


»s2 


®2lens  ^screen  y2f ilm 


H,  *  1.7  lp/mm  (lp  =  line  pairs  at  object  field),  and 
lens 


"film  "  5  2  lp/nnB  • 


screen 


The  lens  dimensions  on  the  3M  Scotchlite  Brand  7610 


reflective  sheeting  determine  V  .  Per 

screen 


manufacturer's  literature,  these  lenses  are  2/1000  inch  in 


diameter.  They  then  constitute  a  resolution  element  equal 

.-2 


to  2/1000  in.  x  25. A  mm/in.  *  5.08  x  10  mm,  or  9.8 
lp/mm. 


Then, 


[_1  • '  5.2*  9.8  J 


-1/2 


=*  1.59  lp/mm  . 


This  number  (H  )  does  not  take  into  account  second-order  effects  such  as 

8 

screen  nonuniformities  and  film  developing  and  is  therefore  an  optimistic 
figure  of  resolution.  The  fact  that  this  number  is  relatively  close  to  the 
empirically  derived  0.8A  lp/mm  lends  additional  credibility  to  the  data. 


*Por  a  7-in.  Aero  Ektar  lens  at  an  f-stop  of  2.5  Optical  Image  Evaluation, 
Proceedings  of  1IBS  Semicontinental  Symposium,  Circular  526,  p.  193. 


DISTANCE  (mm) 


Figure  9.  Impulse  Response  at  Center  Field 

Another  common  reference  figure,  the  minimum  aperture  width,  can  be 
determined  from  the  impulse  response  of  the  system  (Figure  9).  The  aperture 
width  is  the  mimimum  resolvable  spot  diameter  and  is  equivalent  to  the  width 
between  the  half-power  points  of  the  impulse  response.  This  point  is  at  a 
relative  amplitude  of  1.2  on  the  impulse  response  plot  and  corresponds  to  an 
aperture  width  of  0.7  mm  (the  width  of  the  pulse  at  an  amplitude  of  1.2). 

Although  the  limiting  resolution  and  aperture  width  refer  to  absolute 
limits  on  system  performance,  a  considerable  amount  of  research.  Reference  A, 
has  been  done  on  the  required  resolution  for  designated  tasks.  This  research 
is  interesting  because  the  system  end  result  is  an  image  that  must  be  analyzed 
by  an  operator.  Table  II,  which  summarizes  the  results,  indicates  that  the 
system- limiting  resolution  should  be  decreased  (divided)  by  the  factor  of 
Table  II  for  a  given  task.  For  example,  to  specify  the  orientation  of  an 
object,  the  operator  would  require  l.A  line  pairs  per  minimum  target 
dimension.  This  dimension  essentially  reduces  the  limiting  resolution  to  0.8A 
divided  by  l.A  or  0.6  lines  per  millimeter.  The  minimum  target  dimension  for 
orientation  specification  is  the  inverse  of  this  result  or  1.67  mm  (with 
respect  to  the  object  field). 


Using  Royal  X-Pan  (Eastman)  film,  we  obtained  adequate  exposures.  How¬ 
ever,  neither  the  aeroballistics  facility  nor  the  light-source  manufacturer 
had  detailed  light-source  intensity  specifications.  To  determine  the  useful¬ 
ness  of  the  available  light  for  our  purposes  and  for  other  imaging  sensors, 
such  as  television  cameras,  we  needed  to  evaluate  intensity  and  spectral  dis¬ 
tribution;  therefore,  the  light  source  was  brought  to  Los  Alamos  for  that 
evaluation. 

In  the  following  paragraphs,  we  establish  the  measurements  required  to 
evaluate  the  light  source  and  describe  data  collection  and  analysis  and  inter¬ 
pretation  of  that  data. 

a.  Determination  of  Measurement  Requirements.  Radiometric  units  will  be 

used.  In  Figure  10(b),  the  radiant  intensity  1^  of  the  source,  in  watts  per 

steradian,  is  related  to  the  radiant  flux  $  falling  on  an  area  A  of  the 

2  8 

shadow  screen  and  the  irradiance  B  (W/m  )  of  that  area  by 

e 

4>  3  i  Q  a  E  A  watts  ,  (5 

e  e  A  e 

where  is  the  solid  angle  subtended  by  A  in  steradians.  This  relation 

assumes  that  E  is  uniform  over  A.  Then 

0 


A_ 

d2 


(6 


whence 


(7 


The  light  E.  received  at  the  camera  focal  plane  is 
fp 


BeB(9)T  IeR(9)T  2 

Bfp  *  4F2(m  +  l)2  *  4d2F2(m  +  l)2 


In  Equation  (8), 

R(9)  is  the  reflectance  of  the  shadow  screen  at  angle  9  between 
camera  and  light  source; 


(8 


T 


is  the  lens  transmission; 


F 


is  the  aperture  of  the  lens  [that  is,  F  =  f/D, 
where 


f  is  the  lens  focal  length  and 

D  is  the  lens  diameter  (really  the  diameter  of  its  entrance 
pupil)];  and 

m  is  the  magnification  from  object  0  (shadow)  to  image  I  at  focal 
plane  (that  is,  m  =  I/O). 

b.  Parameters  at  Eglin.  The  magnification,  F-stop,  and  transmission  of 
the  lens  at  Eglin  are 

m  =  0.0519,  F  =  2.5,  and  T  =  0.8  . 


To  be  consistent  with  the  values  for  m  and  A  and  to  have  a  18.4-cm 
(7.5-inches)  focal  length,  d  can  be  shown  to  be  3.52  m  (139.2  inches). 

At  least  one  lens  element  will  be  of  lime  glass  for  which  the  transmis¬ 
sion  drops  rapidly  outside  the  360-  to  4000-nm  range. 
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DIVERGENCE  ANGLE  IN  DEGREES 
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(a)  (b) 

Figure  11.  Shadow  Screen  Reflection,  (a)  Directional  reflectance  versus 
divergence  angle  at  0  degree  incidence  of  3M  Ho.  7610  high-intensity 
reflective  sheeting,  (b)  Directional  reflectance  versus  incidence  angle 
at  0.5  degree  divergence  of  3M  No.  7610  high- intensity  reflective 
sheeting. 


c.  Shadow  Screen  Reflectance.  The  shadow  screen  at  Eglin  is  fabricated 
of  3M  Ho.  7610  high- intensity  reflective  sheeting.  This  material  has  a 
reflection  gain  over  standard  diffuse  (Lambertian)  white.  Figure  11(a)  shows 
the  gain  as  a  function  of  divergence  angle  9.  The  reflectance  changes  very 
little  for  oblique  incidence  out  to  20  degrees,  as  Figure  11(b)  shows. 

To  achieve  high  reflectance  gain,  9  should  be  small;  that  is,  the 
camera  and  light  source  should  be  as  close  together  as  possible.  At  Bglin, 

light  source  and  camera  sizes  set  the  minimum  space  between  the  two  at  about 

* 

18  cm.  For  this  space  and  with  3.52  m  for  d,  9  ranges  from  a  minimum  of 
1.2  degree  to  a  maximum  of  3  degrees  for  which  20  <  R(9)  <  100  (if  we 
presume  to  safely  extrapolate  the  curve  of  Figure  11(a)  out  to  3  degrees). 

d.  Focal  Plane/Source  Intensity  Relation  at  Bglin.  Assuming  R(9) 
averages  60  and  the  other  values  determined  in  Equation  (8),  we  have 

E,  =  0.1391  (at  Bglin)  .  (9) 

tp  e 

e.  Source  Intensity  Determination.  A  spectroradiometer  having  well- 
defined  input  geometry  so  that  spectral  intensity  could  be  determined  in  a 
single  set  of  measurements  was  not  available.  Therefore  a  relative  spectral 
distribution  was  determined  using  a  Tracor  Horthem  TH-1710-25  Grating 
Spectroradiometer  equipped  with  a  TN  123-32  head.  Total  energy  output  in  a 
known  geometry  with  a  photodiode  of  known  spectral  sensitivity  was  determined 
with  a  Photodyne  1500XP.  Methods  used  to  extract  the  required  data  are 
described  below. 

(1)  Relative  Spectral  Distribution.  The  light  source  was  perpendicular 
to  the  entrance  slit  of  the  spectroradiometer  and  far  enough  away  that  its 
output  did  not  saturate  the  spectroradiometer  sensor.  The  sensor,  a  linear 
charge-coupled  device  (CCD)  array  located  at  the  output  of  the  grating,  stores 
charges  proportional  to  the  spectral  distribution  for  read-out  after  the  light 
pulse  has  occurred.  The  spectrum  obtained  is  therefore  the  average  spectral 
distribution  during  the  pulse. 


^Camera  and  light  source  cannot  be  made  coincident,  for  then  the  camera  would 
see  only  the  projectile  and  not  the  projectile's  shadow. 
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The  relative  spectrum  S(X)  is  shown  in  Figure  12.  The  ordinate  is 

GI  (X)R(X) ,  where  G  is  an  unknown  calibration  factor,  I  (X)  is  in 
e  _  e 

watts/steradian-nanometers ,  and  R(X)  is  the  responsivity  of  the  CCD  in 
amperes/watt .  The  latter  information  is  available  from  the  manufacturer  in 
terms  of  CCD  quantum  efficiency.  The  responsivity  and  quantum  efficiency 
n(X)  have  the  following  relation  (Reference  3,  page  149). 

R(X),!!XaiXl  (10) 

nc 

The  output  now  becomes 

S(X)  *  G*Xn(X)I  (X)  ,  (11) 

e 

where  G*  includes  the  e/hc  factor.  Values  n(X)  for  the  CCD  are  plotted 
versus  wavelength  in  Figure  13. 


WAVELENGTH  (nm) 


Figure  12.  Relative  Spectra 


WAVELENGTH  (nm) 


Figure  13.  Quantum  Efficiency  of  the  Spectroradiometer 


The  Figure  12  abscissa  was  originally  unsealed.  The  three  dominant  spec¬ 
tral  lines  were  identified  by  comparison  with  known  sources:  770  nm  from 
oxygen,  656  nm  from  hydrogen,  and  501  nm  from  nitrogen.  From  these  value  .e 
scale  was  determined. 

Normalizing  the  responsivity-corrected  relative  spectral  values  *(X)  to 
those  at  a  selected  wavelength  N,  we  get 


S(X)  G'Xn(X)Ia(X)  (12 

S(N)  "  S(N) 


from  which 


!«<>■> 


1 

Xn(X)  x 


S£X1 

S(M> 


SiN l 
G* 


(13 


The  spectrum,  corrected  for  responsivity  wavelength  dependence,  was  normalized 
for  H  a  710  nm,  where  the  amplitude  equals  71.7.  The  result  is  plotted  in 
Figure  14. 


WAVELENGTH  (nm) 


Figure  14.  Corrected  Normalized  Spectrum 

In  Equation  (13),  all  factors  are  known  except  6* ,  which  shall  be 
determined  from  the  total  radiant  energy  measurement  described  in  the 
following  section. 

(2)  Total  Radiant  Energy.  A  Photodyne  1500  XP  radiometer  measured  total 
radiant  energy.  The  light-source  axis  was  aligned  perpendicular  to  and  cen¬ 
tered  on  the  photodiode  sensitive  surface.  The  sensitive  surface  (diameter  a 
3  2  mm)  was  located  at  a  distance  s  (s  =»  65.4  cm)  from  the  light  source.  The 
radiometer  output  is  wavelength  dependent  and  is  in  P(\)  V/w  if  only  the 
sensitive  area  (or  less)  is  irradiated.  Uhen  more  than  the  sensitive  area  is 

irradiated,  as  in  this  measurement,  the  radiometer  output  units  should  be 

2 

changed  to  a  new  value  P'(X)  in  volts/watt  ta  /4,  where  a  is  now  in 
meters . 
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Figure  IS.  Radiometer  Output 


The  output,  a  time-dependent  voltage  corresponding  to  the  light  pulse,  is 


P’(X)i  (X,t)Q  dX 
e  a 


where  Q  is  the  solid  angle  subtended  by  a  and  i  (X,t)  is  the 
a  ~  e 

time-dependent  spectrum  of  light  received.  Figure  IS  shows  voltage  pulse. 

The  spectroradiometer  measure  I  (X)  is  the  time  integral  of 

6 

i0(X,t).  Therefore,  for  comparison  the  radiometer  output  must  be  inte¬ 
grated.  (This  comparison  eliminates  questions  ~  .at  the  time  dependence  of 
the  spectral  distribution.) 


When  we  integrate.  Equation  (14)  becomes 


/ 


v(t)dt 


P’UHett.tJQa  dXdt 


2 

P((X)Xt<X)Q^&t  dX 


where  Z  (X)  is  the  average  spectral  distribution  of  Equation  (13)  and  At 

e 

is  the  pulse  length.  We  insert  Equation  (13)  in  Equation  (15): 
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v(t)dt 


Q 
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•_1_  Sixi]  [sgjQAtj 
Xn(X)  S(N)  J  L  G*  J 
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where  the  bracketed  value  is  the  normalized  spectrum.  But  P* (X) 
2  2 

and  Q  »  «a  /4s  ;  hence.  Equation  (16)  becomes 
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Figure  16.  Normalized  Radiometer  Sensitivity 


A  normalized  curve  P(X)  referred  to  the  sensitivity  at  810  nm,  where  P(X)  = 
1  mV/pW  (Figure  16).  Then  P(X)  =  P(810)  x  P(X)/P(810),  and 
Equation  (17)  becomes 


/v(t)dt  - 

0 


S(mP(810)At 

32g> 


-PCX) 

P(810) 


I 


SIX).-] 


Xn(X)  S(H)J 


dX 


(18 


In  Equation  (18),  all  values  are  known  except  G* .  The  two  brackets  are  values 
for  the  normalized  radiometer  response  and  the  normalized  spectrum.  When  we 
solve  for  ,G’ , 


*2 

f  PCXi 

r  i 

scxn 

/  P(810) 

LXn(X) 

S(H)J 

dX 


S(M)P(810)dt 


/ 


At 

v(t)dt 
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The  value  of  the  wavelength  integral  for  X^  =»  432  nm  to  = 

792  nm  is  240.  The  value  of  the  time  integral  is  0.717.  For  the  radiometer 
sensitivity  scale,  we  used  P(810)  =  1  mV/pW  or  103  V/W.  When  we  substi¬ 
tute  s  and  the  other  values  now  determined,  G'  »  16.8. 

Referring  to  Equation  (13),  we  now  can  describe  the  spectrum  I  (X)  as 

e 


Ie(X)  S(M)  x  the  responsivity-corrected  and  normalized  spectrum 
G' 

s  4.27  x  the  normalized  spectrum  of  Figure  10. 


f.  Illumination  Conclusion.  The  spectral  radiant  intensity  of  the  light 
source  is  now  determined.  For  use  with  any  imaging  sensor,  all  that  is  needed 
is  the  responsivity  of  that  sensor.  The  term  Ifl(X)  is  a  general-purpose 
statement  of  light-source  spectral  output  and  will  be  used  further  in  system 
design  and  specifications.  Information  derived  from  the  resolution  analysis 
can  be  coupled  with  this  result  to  encompass  system  parameters  and  maximize 
S/H  ratio  for  optimum  system  performance. 


5. 


El  SYSTEM  CONSIDERATIONS 

A  number  of  possibilities  exist  for  replacing  film  cameras  with  El  hard¬ 
ware;  each  possibility  will  match  imaging  requirements  established  in  previous 
sections.  However,  processing  strategies  and  associated  hardware  will  vary 
dramatically  for  each  method.  In  all  cases,  electronics  must  be  associated 
with  the  El  unit  to  control  its  operation  and  store  the  image  data.  The 
storage  function  is  necessary  to  permanently  archive  the  image.  El  devices 
are  only  capable  of  retaining  image  data  for  a  short  time  before  data  degrada¬ 
tions  occur;  therefore,  associated  electronics  must  be  able  to  store  the  image 
data  at  video  or  near-video  rates.  Figure  17  shows  a  configuration  for  an  El 
subsystem  that  would  replace  the  present  film  camera.  The  external  communica¬ 
tion  link  provides  a  data  path  from  each  El  subsystem  to  a  central  data  col¬ 
lection  computer  (host  computer)  as  shown  in  Figure  17.  The  host  computer  is 
responsible  for  collecting,  maintaining,  and  analyzing  the  images  from  each  El 
subsystem.  For  each  image,  the  analysis  would  consist  of  finding  the  critical 
points  and  the  projectile  and  establishing  the  coordinates  of  those  points. 

The  coordinate  would  be  in  respect  to  the  reference  beads  at  the  top  and 
bottom  of  an  image.  In  addition,  the  host  computer  would  control  each  El 
subsystem  for  pre-  and  post-shot  operations. 

a.  El  Subsystem.  The  purpose  of  the  El  subsystem  is  to  record  an 
image.  The  operator  normally  has  access  to  the  images  after  the  host  computer 
has  read  the  images  from  each  subsystem  (or  station).  The  major  components 
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Figure  17.  The  El  Subsystem  Acquires  and  Temporarily  Stores  the  Image 
Data.  The  host  computer  collects  all  images  over  a  high¬ 
speed  data  link  for  presentation  to  the  operator. 


(Figure  17)  are  the  El  device  (or  camera),  control  logic,  and  data  buffer. 
During  a  typical  test  shot  control  sequence,  control  logic  would  accept  com¬ 
mands  from  the  host  to  prepare  the  substation  (that  is,  the  camera  scrub 
cycle,  etc.).  The  existing  infrared  detector  system  would  activate  the 
camera.  Immediately  after  the  projectile  is  imaged,  the  control  logic  trans¬ 
fers  the  image  from  the  camera  to  the  data  buffer.  The  data  remain  there 
until  the  host  accesses  the  subsystem. 

Several  options  exist  at  the  subsystem  level  and  several  types  of  cameras 
can  be  used.  (They  are  discussed  in  paragraph  F)  In  each  case,  the  rate  at 
which  data  can  be  read  is  variable.  To  maintain  resolution,  data  must  be 
accessed  as  quickly  as  possible.  Currently,  available  control  and  data  buffer 
hardware  will  operate  at  approximately  13  million  bytes  per  second.  This 
means  that  an  image  will  be  read  from  the  camera  within  approximately  one-half 
second,  which  is  acceptable  to  slow-scan  camera  technology. 

The  control  logic  and  data  buffer  can  be  implemented  with  standard,  com¬ 
mercially  available  products;  custom-built  hardware;  or  some  combination  of 
both.  Using  commercial  products  is  expensive,  but  such  use  minimizes  design 
risks,  and  the  products  are  usually  backed  by  warranty  and  a  maintenance 
organization.  Custom  hardware  can  be  much  cheaper  when  purchased  in  quantity 
and  can  perform  better,  but  usually  it  is  not  backed  by  a  maintenance  organi¬ 
zation  (that  is,  users  must  do  their  own  maintenance). 

For  the  aeroballistics  range,  custom  hardware  would  be  much  cheaper  than 
a  commercial  product  because  the  hardware  would  be  designed  specifically  for 
the  El  subsystem  and  therefore  would  require  less  hardware  than  a  general- 
purpose  commercial  unit;  and,  Los  Alamos  is  not  profit  structured  and  there¬ 
fore  the  hardware  could  be  bought  at  cost. 

Primarily,  hardware  optimization  would  be  developed  around  data  compres¬ 
sion  techniques.  These  methods  will  work  expecially  well  in  the  BI  subsystem 
because  an  image  is  composed  of  a  projectile  and  two  Kevlar  lines;  all  remain¬ 
ing  data  could  be  suppressed,  thereby  drastically  reducing  the  size  of  the 
data  buffer.  The  compression  techniques  used  would  be  the  so-called  error 
free  methods  so  that  the  image  resolution  would  be  unchanged. 

b.  Host  Processor.  The  host  processor  will  act  as  a  data  collection 
point  for  all  of  the  El  subsystems.  The  processor  should  be  a  high-speed, 
general-purpose  micro-  or  minicomputer  with  a  large  mass-storage  device  for 
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storing  images  and  a  tape  drive  for  permanent  data  archival.  A  high- 
resolution  graphics  terminal  would  be  needed  to  display  the  images,  and  an 
image-quality  hardcopy  unit  may  be  desirable  to  produce  prints  of  the  images. 

The  host  processor  could  work  in  a  number  of  modes.  The  simplest  mode 
would  interface  with  existing  equipment  to  get  preshot  timing  information  to 
ready  the  EX  subsystems.  In  this  case,  the  host  could  interface  with  the  ASIA 
controller  so  that  control  currently  associated  with  opening  the  camera  shut¬ 
ters  c  ,uld  be  used  to  initiate  commands  from  the  host  to  the  subsystems  to 
prepare  for  the  shot.  A  more  complicated  route  could  consist  of  an  integrated 
approach  to  the  aeroballistics  facility.  This  mode  would  take  advantage  of 
some  of  the  power  of  a  general-purpose  computer  to  collect  not  only  image  data 
but  also  the  data  now  collected  by  HKRME  (including  timing  data).  In 
addition,  control  functions  now  done  by  ASIA  could  be  done  by  the  host. 
Therefore,  with  this  approach,  the  same  computer  that  performs  the  image 
processing  could  replace  the  current  HP  computer  and  the  ASIA.  This 
integrated  approach  would  ease  maintenance  requirements  and  provide  a 
simplified  operator  interface,  but  it  would  require  more  design  engineering 
than  the  simpler  approach. 

The  operator  will  command  all  of  the  EX  system  functions  from  the  host 
console.  After  a  shot,  the  host  will  access  the  BI  subsystems  in  a  sequential 
fashion.  Image  data  will  be  transferred  over  a  very  high  speed  data  link  and 
be  stored  on  disk.  The  time  required  to  transfer  100  images  over  the  data 

link  depends  on  the  speed  of  the  link  and  whether  data  compression  is  used. 

9 

We  recommend  that  data  compression  be  used;  otherwise,  approximately  10 
bytes  of  data  must  be  transferred  and  stored  on  the  host.  After  transfer  is 
complete,  the  operator  may  start  analyzing  the  images.  We  envision  that  the 
operator,  using  a  light  pen,  will  be  able  to  roughly  indicate  critical  points 
on  the  image,  as  shown  in  Figure  18.  The  host  will  use  these  approximate 
indications  to  find,  by  means  of  advanced  pattern  classification  techniques, 
the  exact  coordinates  of  the  critical  points.  In  other  words,  the  operator 
must  indicate  the  approximate  locations  of  the  critical  points  within  a 
certain  error  radius,  as  shown  in  Figure  19.  The  host  could  then  display  the 
numerical  value  of  the  exact  point  locations  along  with  cross-hair  indications 
of  where  the  final  points  are  on  the  projectile.  This  feature  would  reduce 
the  time  required  to  process  each  image,  but  more  importantly,  it  will  give 
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the  resulting  set  of  data  from  all  stations  a  consistency  that  cannot  be 
achieved  through  manual  operations.  As  a  result  of  this  consistency,  the 
effective  accuracy  of  the  system  is  increased. 

When  the  operator  has  processed  all  images,  the  data  can  be  put  on  tape 
in  a  format  compatible  with  the  existing  method  and  then  transported  to  a 
larger  computer  for  further  analysis  as  is  presently  done.  The  image  data 
could  also  be  put  on  tape  and/or  printed  on  an  image-quality  hardcopy  unit. 
Quick-look  data  could  also  be  made  available,  such  as  a  movie  mode.  In  this 
mode  the  sequence  of  100  images  would  be  aligned  and  displayed  in  rapid 
sequence,  thus  giving  the  researcher  a  quick  idea  of  the  projectile  flight 
path. 

c.  Maintenance.  Commercial  products  such  as  the  host  computer  can  be 
maintained  either  by  company  service  representatives  or  by  ARF  personnel. 
Routine  maintenance  like  tape  head  cleaning  and  camera  adjustments  would  be 
reasonable  tasks  for  ARF  personnel.  In  fact,  it  seems  reasonable  that  the 
current  level  of  personnel  at  the  ARF  could  deal  with  a  majority  of  the 
problems  that  may  arise.  Modern  computer  systems  typically  have  on-board 
diagnostics  that  take  some  of  the  guesswork  out  of  troubleshooting.  In  addi¬ 
tion,  a  certain  amount  of  maintainability  can  be  built  into  the  BI  system. 

For  example.  Figure  20  demonstrates  how  local  BI  subsystem  adjustments  can  be 
made.  A  technician  may  do  maintenance  by  connecting  a  portable  monitor  to  a 
particular  BI  subsystem.  The  monitor  would  display  exactly  what  the  camera 
sees  so  that  adjustments  in  the  focus  or  camera  alignment  could  be  made. 

6.  TELEVISION  REQUIREMENTS 

In  television  systems  the  optical  system  forms  the  image  on  a  photoactive 
surface  (photocathode),  where  electrical  charge  released  is  read  out  by  a 
scanning  electronic  beam  in  all  tube-type  image  sensors  (vidicons,  etc.)  or  by 
a  charge  transfer  process  in  photosensitive  array  devices  [CCOs  and 
charge- injected  devices  (CIDs)].  Initially,  our  discussion  will  be  devoted  to 
the  tube-type  image  sensor;  the  terminology  developed  will  be  directly 
applicable  to  the  imaging  array  devices. 

The  system  of  scanning  lines  is  known  as  the  raster.  Except  for  special 
uses  that  are  not  applicable  here,  the  scanning  lines  are  straight  and  paral¬ 
lel.  Interlaced  scanning,  employed  in  the  television  industry  to  reduce 
flicker,  will  not  be  employed  here.  Adjacent  lines  will  be  adjacent  in  time. 
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The  resolution  perpendicular  to  the  scanning  lines  is  determined  by  their 
spacing,  the  reading  beam  profile,  and  the  point-spread  function  of  the  optics. 
The  resolution  along  a  scanning  line  determines  the  required  video  bandwidth. 

a.  Raster  Design.  The  raster  will  read  a  rectangular  area  on  the  image 
tube  that  corresponds  to  the  field  of  view  at  the  object  plane.  The  scanning 
lines  can  be  along  either  dimension  of  the  raster.  Both  schemes  will  be 
examined  so  that  the  one  requiring  the  lower  video  bandwidth  (and  therefore 
less  video  noise)  may  be  selected.  The  actual  raster  dimensions  at  the  image 
tube  will  be  determined  by  the  photocathode  dimensions,  which  in  turn  will 
determine  the  lens  to  be  employed.  The  image  tube  cannot  be  selected  until 
the  aspect  ratio  and  vertical  and  horizontal  resolution  requirements  are 
determined.  Therefore,  these  measurements  will  be  determined  in  the  object 
field. 

(1)  Object  Field  Dimension.  The  object  field  at  the  reflecting  screen 
is  somewhat  smaller  than  the  screen  and  is  determined  by  the  exposed  part  of 
the  film  in  the  camera.  The  relation  between  the  two  sets  of  dimensions  will 
be  determined  by  the  camera  magnification  m.  The  distance  between  the  fiducial 
beads  at  the  screen,  91. 44  cm,  was  used  as  a  reference.  On  one  test  film 
employed  to  determine  the  HTF  and  the  line-spread  function,  the  images  of  these 


Figure  18.  Analysis  of  Imaged  Data  at  the  Host  Computer  Console. 


Pigure  19.  Tolerance  Radius  in  Critical  Point  Selection 

beads  were  spaced  4.746  cm  apart.  The  image-to-object  (I/O)  magnification, 
that  is,  m  =  I/O,  was  0.0519.  (This  magnification  varied  somewhat  from  camera 
to  camera.) 

The  horizontal  and  vertical  dimensions  of  the  exposed  part  of  the  film 
were  5.842  cm  and  12.09  cm,  respectively.  Hence,  in  the  object  field, 

H  3  5.842/m  3  112.6  cm  or  1126  mm,  and 
V  »  12.09/m  -  232.9  cm  or  2329  mm. 


Figure  20.  Local  Maintenance  on  an  El  Subsystem. 


Figure  21.  System  Aperture  Width. 


The  aspect  ratio  V/H  will  be  2.069:1  if  the  scanning  line  is  along  the 
horizontal  dimension  or  0.483:1  if  it  is  along  the  vertical  dimension. 

(2)  Object  Field  Resolution.  The  point-spread  function  is  the  image 
formed  by  the  system  from  a  point  source  of  light.  As  two-point  sources  of 
light  are  brought  closer  together,  the  limit  of  the  system's  ability  to 
resolve  the  two  is  reached  when  their  spacing  Ax  is  equal  to  the  half -power 
width  of  the  point-spread  function.  As  shorn  in  Figure  21,  the  resulting 
signal  (the  dotted  line)  is  one  television  line.  In  paragraph  II. 3. b. ,  Ax 
was  determined  to  be  0.7  mm  in  the  object  field. 

The  number  of  cycles  required  in  each  direction  will  then  be 


N 


H 


1126 

0.7 


s  1608  television  lines  and 


2329 

0.7 


=  3329  television  lines. 


(3)  Bandwidth  Requirement.  Each  scan  line  will  require  a  time  ts  for 
completion.  As  noted  above,  the  bandwidth  is  determined  by  the  resolution 
along  the  scan  line;  that  is,  the  bandwidths  for  horizontal  and  vertical 
scanning  are 


BW  -  H/ t  and 
H  H  s 


BW  =  H  /t 
V  V  s 
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The  scan  times  may  not  be  the  same,  however.  They  will  be  determined  by  the 
total  time  allowed  to  read  the  raster  and  by  the  time  required  for  interline 
pulses.  Because  only  a  single  frame  (the  raster)  is  to  be  read,  the  normal 
interframe  synchronizing  pulses  will  not  be  required. 

Standard  television  interline  synchronizing  pulses  require  a  time  of 
11.11  ps.  If  we  assume  a  synchronizing  time  of  11  ps,  the  time  per  line 
t^  will  be 

t„  =  t  +  11  x  10~6  . 

1  s 

For  a  total  of  H  lines  in  the  raster,  the  total  time  per  raster  t  will 

K  K 

be  NRt^.  The  resolution  perpendicular  to  the  scan  lines  will  be  less 
than  the  number  of  scan  lines  by  a  factor  known  as  the  Kell  factor. 

Reference  5.  This  factor,  determined  by  the  raster  line  spacing,  the  reading 
beam  profile,  and  the  optical  point-spread  function,  was  found  experimentally 
to  range  from  0.53  to  0.85.  These  factors  included  the  lines  required  for 
interframe  processing  in  standard  television.  To  use  the  modem  stable  raster 
system,  which  does  not  require  interframe  lines,  we  assume  a  conservative 
value  of  0.85.  For  a  required  resolution  of  (perpendicular  to  the 
scanning  lines) , 

H_  =»  11  /0.85 
R  p 

Then  the  total  master  osad-out  time  tR  will  be 

t_  *  H  t. /0.85  . 

R  pi 

If  we  assume  a  raster  read-out  time  of  0.3  s,  the  time  per  line  will  be 

t.  =*  0.255/11  s  . 

I  P 

From  the  above  relations  the  following  table  of  values  can  be  developed 
for  horizontal  and  vertical  scanning. 


Scanninx 

s  (vs) 

Bandwidth 

Horizontal 

3917 

54.1 

29.7  KHz 

Vertical 

1892 

123.8 

26.9  KHz 

Although  the  bandwidth*  differ  little,  the  smaller  bandwidth  of  the 
vertical  scanning  system  will  reduce  the  noise  power  by  about  10  percent. 

b.  Selection  of  Image  Sensor.  For  perfect  geometry  the  ideal  sensor 
would  be  a  photosensitive  array  device  such  as  a  CCD  or  a  CID.  At  the  present 
time,  none  of  the  commercially  available  products  meets  the  resolution 
requirements.  However,  as  part  of  project  Galileo  and  the  space  telescope, 
the  National  Aeronautics  and  Space  Administration  (NASA)  is  funding  work  on 
high-resolution  array-type  devices.  A  summary  of  known  work  on  imaging  arrays 
is  presented  in  paragraph  6.e.,  along  with  some  ideas  for  expanding  the 
effective  resolution  of  present  devices. 

Among  the  tube-type  image  sensors,  only  the  return  beam  vidicon  exceeds 
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the  BI  resolution  requirements.  Reference  6.  This  tube  is  not  available 
commercially.  The  Westinghouse  WX-5156,  designed  for  document  reading,  is 
capable  of  reading  3200  lines/scan;  its  HTF  is  5  percent.  If  we  extrapolate 
to  3  percent,  as  Johnson  advises.  Reference  7,  the  resolution  would  be  3746 
lines/scan.  The  imaging  optic,  light  requirement,  dynamic  range,  and  expected 
noise  for  the  WX-5156  are  described  below. 

c.  Raster  Dimensions.  The  useful  diameter  of  the  WX-5156  image  tube  is 
35.6  mm  (1.4  inch).  A  raster  having  an  aspect  ratio  of  2.069:1  must  be  fitted 
within  this  circle.  The  raster  horizontal  and  vertical  dimensions  are  h  and  v 

2  2  2 
h  +  v  =  35.6  and 

v/h  =  2.069 

from  which 

h  -  V 35.6Z/(2.0692  +  1)  *  15.47  mm  and 


v  »  32  mm 


In  the  process,  extra  vertical  resolution  is  obtained.  The  scan  line,  which 
had  a  3746  line  resolution,  was  for  a  3:4  aspect  ratio,  where  the  scan  dimen¬ 
sion  was  28.44  mm.  Therefore,  the  available  resolution  along  the  scan  line 
will  now  be  32/28.44  x  3746  or  4215  lines. 

Two  questions  must  now  be  answered:  Has  the  extra  resolution  along  a 
scan  line  been  obtained  at  the  expense  of  resolution  perpendicular  to  the  scan 
linet  Does  the  bandwidth  have  to  be  increased  to  take  advantage  of  the  avail¬ 
able  resolution  along  a  scan  line? 

The  scanning  beam  profile  determines  the  maximum  useful  scan-line  density 
and  the  resolution  along  a  scan  line  when  resolution  is  not  bandwidth  lim¬ 
ited.  The  half-width  Ab  of  the  scanning  beam  profile  is 

scan  line  length  32  , 

^  *  lines/scan  3  4215  3  7.59  x  10“  mm  . 

The  maximum  useful  scan  line  density  is 


max 


15.47 

7.59  x  NT3 


2038  . 


This  value  exceeds  the  required  value  for  H  (for  vertical  scanning)  by  8 

tc 

percent . 

If  the  full  resolution  capability  is  to  be  attained,  the  bandwidth  will 
have  to  be  increased;  that  is, 

BW  »  4215/123.8  x  10~6  *  34.0  MHz  . 
r 

However,  note  that  the  increased  resolution  will  be  only  along  the  scan  line 
(vertical  dimension).  The  resolution  perpendicular  to  the  scan  lines  (the 
horizontal  dimension)  will  be  limited  to  the  value  determined  in  Sec.  F.l.b. 

d.  Imaging  Optic.  The  required  focal  length  can  be  determined  from  the 
magnification  required  and  the  lens-tu-object  field  dimension.  From  the  thin 
lens  equation,  we  can  show  that 

«  om 

1  *  m  ’ 
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where  f  is  lens  focal  length,  o  ■  lens-to-object  distance,  and  m  *  the  I/O 
ratio.  From  paragraph  4.b.,  o  =  3520  mm.  The  I/O  ratio  is  32/2329  or 
0.0137.  Then, 


f 


3520  x  0.0137 
1.0137 


47.7  mm  , 


a  value  within  95  percent  of  the  standard  50-mm  focal  length  lens  employed  in 
35-mm  cameras.  For  f  a  50  mm,  m  *  0.0144.  A  small  compensating  adjustment  in 
the  vidicon  deflection  optics  should  be  possible  without  degradation  of 
geometry  or  resolution. 

At  the  focal  plane,  a  resolution  of  4215/32  or  132  television  lines  is 
required  for  full  bandwidth  or  3329/32  =  104  lines  for  minimum  bandwidth.  The 
latter  is  easily  obtained  with  modem-design,  off-the-shelf  lenses.  However, 
the  lens  resolution  must  exceed  the  complete  system  resolution  by  a  sufficient 
margin  so  that  lens  resolution  does  not  degrade  system  resolution  appre¬ 
ciably.  Component  resolution  controls  the  system  resolution  Hg,  as  the 
following  relation  shows: 


Three  known  components  are  the  lens  H. ,  the  screen  11  ,  and  the 

l  sc 

vidicon  H^.  In  the  object  field, 

V  *  9.8  cycles/mm,  a  value  determined  by  the  size  of  the  screen  lens; 
sc 

Nv^  *  mil  a  0.95  cycles/mm,  the  maximum  resolution  available  from  the 
vidicon;  and 

Hg  ■  (104/2)m  »  0.75  cycles/mm,  the  minimum  allowable  system 


resolution. 


Solving  for  11^  gives  the  minimum  admissible  lens  resolution: 


"life  -(~h  *  -H 

L*.  V*  *v  *.Ji 


0.588 

0.0144 

*  40.8  cycles/mm  (with  respect  to  the  focal  plane) 


This  value  will  require  a  high-quality  lens,  which  can  be  purchased  from  a 
number  of  commercial  sources. 

Light  Requirement 

The  minimum  light  will  be  that  amount  producing  a  signal  at  the  limiting 
resolution  spatial  frequency,  which  is  recognizable  in  the  prevailing  video 
noise.  Therefore,  video  noise  must  be  determined  first.  The  final  determina¬ 
tion  must  compare  signal  and  noise  electron  inputs  with  the  video  amplifier 
system — a  problem  in  recognition  statistics  (Reference  3,  Chapter  8). 

Video  Boise 

With  vidicon  image  tubes,  noise  in  the  first  stage  of  the  video  system 
will  dominate.  A  good  video  front  end  will  have  noise  current  Lnput  equiva- 

_9 

lent  to  about  7  x  10  A  for  a  30-HHz  bandwidth.  For  white  noise  this  input 

_9  .  _i2 

is  7  x  10  /TJ0  MHz  or  1.3  x  10  A/Hz.  For  maximum  resolution,  the 
bandwidth  required  is  34.0  MHz.  The  equivalent  video  noise  input  would  be 


1.3  x  10 
7.6  nA  . 


x  \34.0  x  10 


This  video  noise  input  constitutes  a  noise  electron  arrival  rate  of 
A  in  7.6  x  10"9  C/s 


1.6  x  10- 


C/electron 


»  4.74  x  10  electrons/s 


*.  •.  *.  \  A  *,  r 


•-*  -•  ■■  ■•  -  •  -  ’  oW\ 
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Picture  element  Time 

The  time  At,  which  is  required  to  read  out  each  picture  element  at 
maximum  resolution,  v'U  be 


fcs  123.8  x  10-6  „  ..-9 

At  *  r  * - mi -  29  x  10  8  * 

v 

During  this  time  the  noise  electrons  collected  will  be 


n„  =  rifjAt  =  A. 74  x  1010  x  29  x  10-9  =  1392  and 


4nn  *  "n  *  37 


The  variation  in  the  noise  electrons  will  be  less  than  3  percent;  therefore 
average  values  will  be  employed.  Recognition  statistical  treatment  is  not 
required. 

S/N  Ratio  at  3%  MTF  and  Required  Light  Level 


Por  standard  broadcast  television,  the  S/N  ratio  at  the  limiting  resolu¬ 
tion  can  be  less  than  one  because  the  eye  integrates  the  displayed  scene  for 
about  0.2  s.  This  is  time  to  correlate  out  the  noise  over  6  frames  (or  12 
fields).  In  this  application  a  single  frame  will  be  read.  Visual  experience 
while  looking  at  an  A  scope  oscilloscope  display  of  a  single  line  of  video 
shows  that  a  signal  having  S/N  *  1  cannot  be  seen,  whereas  for  S/N  *  3  the 
signal  can  be  seen.  At  S/N  =>  2  the  signal  can  be  seen  occasionally.  There¬ 
fore,  let  S/N  .  ■  3  (that  is,  n  x  MTF  (0.03)  >  3n  ],  where  n  are 
min  s  ns 

the  signal  electrons  that  would  be  collected  during  a  picture  element  time 

At  if  the  MTF  was  1.  Then  n  *  3n  /0.03  =*  lOOn  ,  providing  a  dynamic 

o  n  ^ 

range  of  at  least  100:1  of  which  the  system  requires  ng  a  1.39  x  10 ‘ 
electrons. 

Selection  of  Photoconductor 

Four  different  photoconductors  are  available  with  the  WX-5156  (Table  3). 
They  differ  primarily  in  spectral  response  and  image  degradation  with  time 
after  exposure.  All  show  a  small  rise  in  output  until  about  2  s  after 
exposure.  (This  increase  is  probably  due  to  the  so-called  soaking  phenomena 
caused  by  release  of  charges  from  deep-lying  sites  within  the  photo¬ 
conductor.)  How  much  time  elapses  after  each  exposure  will  depend  on  the 


associated  computer  system  architecture.  Therefore,  tentatively,  the  photo¬ 
conductor  with  minimum  image  degradation  will  be  selected.  This  is  photo¬ 
conductor  B  that  degrades  in  limiting  resolution  by  32  percent  at  100  s  after 
exposure  (Figure  22).  Unfortunately,  this  is  a  less  sensitive  photoconductor 
(Figure  23).  However,  the  spectral  response  of  photoconductor  B  (Figure  24) 
is  trail  matched  considering  the  light  source  and  the  decay  characteristics. 

It  peaks  at  4250  A,  and  the  spark-gap  light  peaks  in  the  blue  spectral 
region  (Figure  14). 

Signal  Electrons  n8 

The  signal  electrons  will  be  collected  during  the  spark-gap  flash. 

Except  for  a  small  increase  as  a  result  of  soaking  and  a  small  loss  between 
exposure  time  and  collection  time,  the  number  of  electrons  will  be  propor¬ 
tional  to  the  total  radiance  received  from  the  light  flash  and  the  spectral 
match  between  the  light  and  the  photoconductive  material  B.  The  spectral 
response  R(\)  is  given  in  amperes/watt.  The  number  of  signal  electrons 
released  per  unit  area  Hg  at  the  photocathode  will  be 


A. 

/ 


R(X)Efp(X) 


dX  signal  electrons/m 


where  Bfp  is  the  radiant  energy  received  at  the  focal  plane  (the  photo¬ 
cathode)  .  is  related  to  the  radiant  intensity  of  source  Ig  by 

Equation  (8),  paragraph  4. a.  Then, 


„  ,  —  g-Wl - -  r2  R(x)le(X)d\  . 

4d2eF2(m  +  l)2  J  e 

In  paragraph  4.e.,  Ig(X)  «  4.27  x  the  normalized  light-source  spectrum  of 
Figure  14.  The  new  value  for  m  is  0.0144.  Using  the  same  values  for  the 
remaining  terms  outside  the  integral  as  those  used  in  paragraphs  4.b.  and 
4.c. ,  we  reduce  the  terms  to  9.4  x  1017.  Then, 


>s 


9.4  x  1017 


R(X)Ie(X)dX  . 
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STORAGE  TIME  (s) 

Figure  22.  Storage  Characteristics  of  Slow-Scan  Vidicons 

Between  4000  and  6000  A,  the  integral  reduces  to  4  x  10  Then, 

16  2 

B  »  3.6  x  10  electrons/m  ,  but  n  ■  H  AA,  where  AA  is  the 
s  s  s 

minimum  picture  element  area.  In  the  15.47-  by  32-mm  raster,  there  are  a 
maximum  of  4215  x  1608  picture  elements.  Therefore, 

15.47  x  32 

AA  =  4215  x  1608  *  7.3  x  10” -*  mm^/picture  element 
or  7.3  x  10-H  m^/picture  element 

and 

ns  =«  VSAA  »  3.6  x  10**  x  7.3  x  10”^  >  2.6  x  10*  electrons/ 
picture  element 

and  the  requirement  that  S/ll  =  n  /n  >  100  is  more  then  met.  The  F-stop 

s  n 

of  the  lens  can  be  increased  for  improved  resolution  until  diffraction  limits 
resolution  to  less  than  that  desired. 
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TABLE  3.  CHARACTERISTICS  OF  SLOW-SCAN  VIDICONS 


WAVELENGTH  (ft) 


Figure  24.  Spectral  Response  of  Slow-Scan  Vidicons 

e.  Possible  Camera  Sources.  Table  4  is  a  fairly  comprehensive  list  of 
current  BX  device  manufacturers,  which  includes  current  products  and  products 
planned  for  the  future.  Vote  that  a  VASA/ Jet  Propulsion  Laboratory  (JPL) 
program  is  presently  developing  an  8000  x  8000  element  CCD  array  for  a  space 
telescope  application.  This  will  be  done  by  coupling  64  1024  x  1024  element 
arrays  by  means  of  fiber  optics.  We  have  had  experience  with  this  technology 
from  developing  a  high-speed  BI  system  for  an  Air  Force  missile  fault  analysis 
program.  We  found  the  fiber  optic  coupling  to  be  very  effective,  but  it  is 
expensive. 


SECTION  rn 

CONCLUSIONS  AND  RECOMMENDATIONS 

Tests  were  done  to  derive  the  parameters  governing  operation  and  resolu¬ 
tion  of  the  ARF  high-speed  camera  system.  In  particular,  the  system  transfer 
function  was  found  and  the  spark-gap  flash  light  output  was  quantized.  In 
addition,  our  findings  were  presented  in  support  of  the  high-speed  electronic 
imaging  concept. 

The  system  transfer  function  includes  the  camera,  camera  optics,  film, 
film  developing,  and  effects  imposed  by  the  light  source,  light  path,  and 
reflective  materials.  The  transfer  function  (or  MTF)  is  an  industry  standard 
for  specifying  such  systems  and  gives  a  concrete  definition  of  the  camera 
system  resolution.  From  the  MTF,  a  3  percent  response  point  has  been 
calculated  as  an  industry  standard  reference  point.  At  this  point,  the 
minimum  object  field  dimension  that  can  be  resolved  is  0.7  mm:  over  the 
object  field  of  view,  the  system  is  capable  of  resolving  a  matrix  of  data 
points  that  is  3329  elements  high  and  1608  elements  wide.  An  El  subsystem 
must  approach  this  resolution. 

Two  types  of  El  devices  approach  this  resolution:  an  FPS  vidicon  and  a 
large-target-area  slow-scan  vidicon.  The  FPS  vidicon  can  exceed  the  required 
resolution,  and  slow-scan  cameras  are  currently  operating  at  maximum  perform¬ 
ance  at  this  resolution.  The  slow-scan  camera  is  more  desirable  because  of 
linearity  considerations. 

Considering  industry  trends  is  important  for  this  application.  Modem 
solid-state  physics  is  causing  a  technological  revolution  in  the  El  field. 
Solid-state  imaging  devices  are  replacing  vacuum  tube  devices  and  are  allowing 
smaller  size  and  improved  dynamic  response  and  linearity.  However,  solid 
state  has  yet  to  better  the  resolution  or  price  of  traditional  devices. 

Because  solid-state  cameras  use  the  same  technology  as  other  solid-state 
devices  (integrated  circuits,  computers,  etc.),  we  expect  that  performance 
(resolution)  will  increase  rapidly  and  price  will  decrease,  and  we  see 
evidence  of  this  trend  already. 

Analysis  of  the  spark-gap  flash  indicates  the  spectrum  is  fairly  flat 
from  400  to  800  nm,  which  is  compatible  with  most  El  devices.  The  total  light 
output  of  the  unit  was  also  calculated.  A  survey  of  the  light  energy  and 
spectral  sensitivity  of  various  El  devices  indicates  that  the  light  output 
would  probably  be  sufficient  for  most  devices. 


TABLE  4.  CAMERA  SOURCES 


Camera 

Phone 

Contact 

Comments 

Slow— 

Scan 

Video 

Sierra 

Scientific 

408-733-6730 

William  Oyle 
(Marketing) 

Paul  Mongers 
(Sales) 

Harvey 

Bartholomew 

(Engineering) 

Could  have  one  within 

4  mo.  with  1608  x 

3339  res.  and  approp. 
aspect  ratio  and  scan 
rate 

Westinghouse 

607-796-3427 

Rolf  Myer 
(Engineering) 
Pious  Hold 
(Sales) 

1608  x  3329  with  1/2-s 
read-out 

EPS 

V  id icon 

Westinghouse 

High-resolution 

2500  x  2500 

Solid 

State 

Fairchild 

415-858-6149 

Allan  King 

488  x  380  available. 

Will  call  back  on 
advanced  developments 

RCA 

717-397-7661 

Parren 

Gottshall 

(Sales) 

403  x  512  available, 

562  x  512  in  2  years 

Thompson-CSF 

201-438-2300 

John  Mu 1 row 
(Sales) 

576  x  384  available 

GB  Syracuse 

315-456-2832 

-2808 

Dave  Hunter 

380  x  488  available, 

1024  x  1024  in  2-3  years 

Texas 

Instruments 

214-995-2011 

Mr.  Hoseck 

800  x  800  available, 

1024  x  1024  now  in  lab 
for  NASA-JPL  (Jim  Tamand 
213-354-7734) 

Grumman 

Aerospace 

Corp. 

714-660-4200 

Stan  Schwartz 
(Manager  of 
Engineering) 

Stackable  linear  array 
on  thin-film  substrate. 

An  El  device  requires  support  electronics  to  control  camera  functions  and 
temporarily  store  the  image  acquired  by  the  El  device.  An  El  subsystem  is 
composed  of  the  device  and  the  support  electronics.  At  the  ARF,  each  camera 
would  be  replaced  by  an  El  subsystem  that  would  be  linked  to  a  host  computer 
by  means  of  a  single  high-speed  data  link.  The  host  computer  collects  images 
from  the  subsystems  for  presentation  to  the  operator.  Operator  interaction 
coupled  with  advanced  pattern  classification  techniques  is  suggested  to  locate 
critical  projectile  points  with  a  consistent  accuracy. 

Commercial  products  can  be  used  for  the  suggested  system,  but  software 
would  need  to  be  generated  to  provide  the  necessary  functions.  Options  exist 
within  the  El  subsystem  for  custom-built  support  electronics  that  may  greatly 
reduce  the  costs  of  upgrading  the  entire  aeroballistics  range.  The  disadvan¬ 
tage  of  custom  electronics  is  that  usually  they  cannot  be  maintained  by  a 
commercial  service  organization. 

With  this  report  concluding  the  feasibility  study,  it  becomes  clear  that 
BI  can  be  used  for  high-speed  imaging  at  the  aeroballistics  range.  Economic 
feasibility  is  a  question  that  the  funding  agency  must  address.  Because  of 
the  current  technological  climate,  the  economics  of  the  project  will  improve. 
This  advantage  must  be  weighed  against  the  ever  increasing  cost  of  maintaining 
the  current  facility  (caused,  for  example,  by  increasing  film  costs,  downtime 
due  to  component  failures,  and  unavailability  of  parts). 

The  logical  continuation  of  this  project  is  to  build  a  prototype,  a 
state-of-the-art  system  for  high-speed  imaging.  Building  a  prototype  will 
allow  us  to  solve  problems  that  may  occur  in  a  state-of-the-art  design.  Thus, 
El  technology,  as  applied  to  high-speed  imaging,  will  be  in  a  state  of  readi¬ 
ness  as  applications  emerge. 
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